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1 Introduction 

Bogolubov's forest formula [l]-[5] is one of the genistoncs of Quantum Field Theory (QFT), used to handle 
iterative renormalization procedure. In [6] and especially in [7] it was shown that forest formula is actually 
a pure algebraic construction, devised to solve iteratively a broad set of non-linear equation and, remarkably, 
one and the same formula describes everything: from solution to an ordinary quadratic equation to Lamb shift 
and renormalization of the Standard Model. The general presentation of ref.[7] can be a little too abstract and 
does not contain examples. Some examples can be found in [8]-[14], but they arc all at the level of quantum 
field theory and therefore obscure the simple algebraic and tcnsorial structure. A pedagogical review [15] deals 
with adequately simple examples, but it does not pay enough attention to algebra and forest formula. It is a 
purpose of this paper to provide the needed illustrations to the formalism of [7], and this is a natural addendum 
to generic presentation of non-linear algebra in ref.[16]. 

1.1 The problem 

Forest formula provides iterative solution to the following problem. 

• Let F{T) be a given function of some variables T, called times or coupling constants. In QFT applications 
F{T) is usually an effective action (logarithm of partition function) of a physical theory, considered as a function 
of bare couplings, but this interpretation is inessential for the forest formula per se. 

• Let V- be a linear projector onto the linear sub-space of unwanted functions. In renormalization-theory 
context the function F{T) depends also on a cut-off A, and unwanted are functions, which are singular in the 
limit A — > DO. However, again, concrete interpretation of projector V- is inessential for the forest formula. 

• Given F{T) and V- we are ready to formulate the problem. Can one make a change of arguments T in 
F{T), without changing the shape of this function, to put it into the sub-space of wanted functions? In other 
words, the problem is to find a change of variables T ^ T = T + Q{T), such that 

P-{f(t + Q{T)^} =0 (1) 

To make solution of this problem unambiguous, one can impose additional constraint that the counter-term 
Q{T) is pure singular, i.e. lies entirely in the space of unwanted functions: 

r+{Q{T)}=0, (2) 

where V+ = I — V- is projector onto the space of wanted functions^. 

1.2 Forest formula: solution to the problem in terms of Feynman diagrams 

Forest formula provides an explicit solution to the problem (1) & (2), and the solution is in terms of graph theory! 
The crucial observation [7] is that any formal series in non-negative powers of couplings T can be represented 
as a sum over Feynman diagrams, with T standing at the vertices, and this is the most adequate language for 
tensor-algebra studies [16], sec examples below. This means that both the F{T) and the counter-term Q{T) 
can be expanded into sums over graphs: 

F{T) = Y,mmn q{t) = y.q{t)z{t\t) (3) 

r r 

with basic functions Z{T\T), capturing the topology of the diagram F: Z{T\T) is a product of couplings T at the 
vertices with indices (if any) contracted along the links. We assume that projectors 'P± act on the coefficients 

F{T) and Q{T) rather tliaii on Z{T\T). and this is the case in most applications. 

In fsLct, this conditions is somewhat more complicated: -p+|F(r/r)Q(r)| = 0, see sections 2.4.1 and 2.4.2 
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Now iterative solution to the problem (1) & (2) is given by the forest formula: 



F{T/T)Q{T) = -V- { 



F{T)+ Yl ^^(r/7i.--7fe)0(7i)---Q(7fe)[ (4) 

{7lU...U7fe} 



where {71 U . . . U 7^} are all possible box-subgraphs in graph T, -fi G BT, i.e. parts of F lying in a collection of 
non-intersecting "boxes", and r/71 ... 7^ is obtained by contracting all boxes to points. The factor F/F is just a 
single vertex with the same external lines as original F, and Q(F) ^ and is given by eq.(4) only if F{T/T) ^ 0, 
i.e. depending on the shape of the function F{T). The meaning of this criterium is simple: it is the coupling at 
this vertex F{r/r) that is shifted by Q{T). We ascribe a special name "vertex criterium" to this rule, because 
of its significance for our further considerations. 

We refer to [6, 7] for discussion of algebraic theory behind eq.(4), its relation to the group of diffeomorphisms 
of the moduli space of theories and to its dual Connes-Kreimer Hopf algebra on graphs. Instead in what follows 
we consider various examples and applications of eq.(4). 

1.3 The second- level forest formula 

In addition to the first-level forest formula (4), there is a next-level Bogolubov forest formula [l]-[3], involving 
sums over embedded box-subgraphs, sec cq.(9.7) of [7]. It gives directly the answer for the renormalized eflfective 
action F{f) = F{T + Q{T)) = Ffl(F)Z(F|T): 



of r 



The sum is taken over all possible forests Tv-, which are unions of trees for every connected component of F. 
For some connected component of F, Fi, the tree T is build in the following way: 

i) Choose concrete sequence of embedded box-subgraphs {7n}: 

70 = Ti, 71 e BF, 72 e B71 C BF, . . . ,7„ e B7n-i C . . . C B71 C BFi (6) 

where graphs 7^ are not obligatory connected and can consist of s(fc) connected components 7^, . . . , 7^^'^''- But 
for every 7^ should exist such 7^^^, that 7^ G ^Ik-i (condition (6) is understood exactly in this way). 

ii) Making use of this sequence of embedded box-subgraphs, or, equivalently, of the collection of non- 
intersecting boxes, build a decorated rooted tree T. With lower site of each box one associates a vertex of the 
tree. Two vertices are connected by a link, if one of the corresponding boxes lies immediately inside another 
(i.e. there are no boxes in between the two). The root link ends at a vertex, associated with 70 = F. In terms of 
graphs, now a vertex of the tree is associated with every 7^, and there is exactly one link which goes downwards 

(towards the root, i.e. associated with the neighboring bigger box) and connects 7^ to 7^*1']^, and unrestricted 

number of links, going upwards and connecting 7/c to some collection 7^_|_i, . . . , 7^+1*^ C ^k+\- 
The tree product 

5 

■vertices 

of r 



associates with every vertex 7^ an operator {F{\cikhk\)^ {~'^-)^ {ikl ilk+i ' • ■ • ' Ik+ii) ' while projector 
V- acts upwards along the branches of the tree. The root vertex 70 (i.e. a connected component of F) contributes 
just -^^70/(71 • • • • • 7i^°^)) • Arrow over the product sign means that the product is ordered along the branches. 



The problem (1) & (2) implies that function F{T) and projector V- are already given. However, the diagram 

technique still needs to be introduced. Moreover, for given F{T) and "P- this can be done in different ways. It 
deserves mentioning also that iterative nature of forest formula does not mean that there is a small expansion 
parameter: moreover, in renormalization-theory context the relevant expansion parameter is rather large than 
small. 
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Part I 

Algebraic examples 

2 Quadratic equation. A valence-2 vertex and tree diagrams 

Our starting example is 

f F{T) = T + eT'^, 

[ V-{f{e)} = f{e)-f{0), 

i.e. "unwanted" are functions, which depend on e and the goal is to eliminate e-dependence by e-dependent 
shift of variables T — > T + Q^{T), so that F{T) = T + eT^ no longer depends on e. Imposing the constraint (2), 
we reduce the problem to 

f + ef^ = T (8) 

and therefore 

^ ^ Vl + 4eT-l y H^2^(2fc-1)!! ^^^^+1 ^T-eT^ + 2e^T^ - be^T^ + Me^T^ - . . . (9) 

2e ^ (fc + i)! 

Now we are going to reproduce the same series from the forest formula (4). To do this, we need first to 
choose a particular diagram technique. Our first choice will be a single valence-2 vertex T = Tij , where indices 
i,j take a single value i,j = 1 and will be omitted at the beginning. The only possible connected Feynman 
diagrams are chains and circles, both equal to T", where n is the number of vertices. Only two diagrams, with 
n = 1 and n = 2 contribute to F{T). 

2.1 Chain diagrams 

Assume first, that these are two chain diagrams:"^ 

F(^) = l, F(— )=e (10) 

Accordingly F{r /T) ^ only for chain diagrams F. By vertex criterium this implies that the only non-vanishing 
are contributions to Q of chain-diagrams. For thorn the forest formula (4) gives; 

Q[l]^Q(^) = -r.{F{^)}=0, 
Q[2] = Qi~) = -V-{f(~)} = -e, 
Q[3] ^ Q(— ) = + 2F(~)Q(~)} = 

= -2P_{f(^^)Q(^^)} = +2e2 

In the last formula ■F( » « » ) = because of our choice of F{T), and coefficient 2 in from of the next term 
appears because there are two different possibilities to choose a box-subgraph » » inside » » » : | r~»~| » and 



Similarly, three box-subgraphs 

h » » | » - , [«""»} and » | » » » h contribute to the next Q: 

Q[4] = Q( » . ■ . ) = -V- ) (2Q(— ) + = -5e^ 

Note that box-subgraphs » r»~~»l » ' and | r~»~| » » - do not contribute, because contraction of boxes produces 
i/iree- vertex diagrams and such contributions would be multiplied by f ( * * * ), which vanishes for our choice 
of F{T). For continuation of the calculation let us denote the AT- vertex chain diagram through [N]. Then 

0[5] = -P_ [f[2] (2Q[4] + 20[2]0[3]) } = 14e^ 

*In what follows we omit hats over Q(T) and F(r), assuming that they will not be confused with the corresponding functions 
Q{T) and F{T) of T-variables. 
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and, in general, eq.(4) implies 



N-2 



Q[N] = -V- I F[2] 5n,2 + 2Q[N - 1] + ^ Q[k]Q[N -k]\\, N>3 



k=2 



For generating function Q{t) = Z^jv=2*^[-^] t^i^ gives: 

Q = -F[2]{t'' + 2tQ + Q^) (11) 

with F[2] = e, i.e. our familiar (9): 

Qit)+t + e{Q + tf = t, Q{t)+t= ^^^*~^ (12) 
2.2 Introducing coefficients 

In above example the basic functions Z{T\T) entered original function F{T) with coefficients 1 and e. If instead 
we write 

F{T) = a — + eb — — = aT + ehT^, i.e. F(— ) = a, F{-^) = eb, (13) 
then the coefficient a will appear at the l.h.s. of eq.(4): 



aQ{^) = -V-{F{^)}=0, 
aQ(„) = _p_{f(„)} = -eb, 
aQ(— ) = -2P_{f(~)Q(~)} 



so that 

T = T + Q{^)T + Q{~)T^ + Q{~^)T^ + ... = T-^+ \ 



and 



efeT^ 2e%'^T^ \ f ebT^ 2e'^b'^T^ 



2.3 Application of the second-level forest formula 

The second-level forest formula (5) produces the renormalized function F{T) directly, we already know what it 
is from eq.(14) and this answer is very simple: 

F{f) = Fr{T) = aT 

This means that in this case the only non- vanishing output of eq.(5) should be Fr{—^ ) = a, while for all other 

graphs Fii{T) = 0. The first statement is immediately reproduced by (5): there is only one terms at the r.h.s. 
of (5). For more complicated F exact cancelation should occur between different contributions to (5). 
Let us demonstrate this for » » , where one can build the following two trees: 

» » and 



We do not consider trees which have box-subgraphs with one vertex inside, like , because 'P-F{—»- ) = 0, 
also we do not draw external box which corresponds to the net graph. Therefore 

Fr{~) = Fi~) + F{-^)—^i~P.)Fi-^) = F(~) - V.F{~) = (15) 

^ (-^) 

Next is » » » with embedded box-subgraphs (only sub-trees that give non- vanishing contributions to (5) are 



shown): |» » | » ,- F~i~| « - and » rT =i| - ,- » | » »| ' • It is quite obvious that the first two and the second two 
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trees give the same contribution, thus we take into account only the first two graphs and multiply the result by 
a factor of 2: 



= 2|F(- 
Thus Fr(-.— 



1 



-^) vanishes, 

)=- 



i-r-)F{- 



F(— ) 



(16) 



{f(— ) - r_F{~)V_F{~)} 



0, 



because F{* * ) is 'P-negative. 

Now it is easy to understand the way in which cancelations occur: all trees for some particular graph can 
be divided into pairs: the only difference between trees in one pair is that one of them has "external" box, i.e. 
in terms of (5), 7o = 7i = Ti. Due to the presence of additional projector {—V-), the two graphs in each pair 
exactly cancel each other. 

In what follows, we do not give examples of application of second level forest formula in cases where it is 
trivial, because there the same argument can be used to proof cancelation all except some peculiar graphs. 



2.4 Generic projector 

As the next generalization we do not fix concrete e-dependence of the coefficients A(e) and B(e) in 

F{T) = A{e)T + Bie)T^ = Aie) ^ + B{e) ~ i.e. F{^) ^ A{e), F{~) ^ B{e) (17) 

(e.g., they can be generic Laurent series), and concrete choice of projector V- (e.g., it can pick up all terms with 
e'^ and k > ko or, instead, all terms with e'^sinhe or anything else). Eq.(4) works in all situations! Indeed, the 
counter-terms are: 

A{e)Q (^) = (^)l = -A_, 



I.e. 



A{e)Q 



= -V- 



F 



+ 20 



Q 



A_ 
~A 



BAl 
~A^ 



Q[3] = Q (— ) = -jP- {F[2] (2Q[2] + 2Q[2]Q[1])} = 



I bI 

A I A 



BAl 



A_ 
~A 



2 / BAj, 
A I ~W 



BAl 
~A^ 



Q[4] = -jV- (F[2]02[2] + 2i^[2]Q[3] + 2F[2]0[3]0[1]) 



A I A^ 



BAl 
~A^ 



4BA^ 



BA+ fBA 



A^ 



A^ 



Therefore 



BAl 
~A^ 



A 



BA+ fBA"] 



A^ \ A^ 



2^3 

~a' 
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and 



F{f) =Af + BT"^ = A+T + 



BA\ (BA\ 



A^ 



BA+ (BA\ 



A^ 



BA+ (BA\ 



A^ 



A^ 



BA+ BA+ (BA\ 
A^ \ A^ \ A2 



BA+ / BA+ ( BA 



A^ 



A^ \ A^ 




B [BA^, 



A^ \ A2 



Clearly, all terms before every power of T gather into 'P-positive expression: 

. BA+ ( BA 



f BA^ 

F{f)= A+T+{^^ \ r2-2 



BA+ (BAl 



A2 



A2 



A2 



A2 



BA\ 
~A^ 



B (BA\ 



A2 V 742 



2^4 _ 

+ 



The r.h.s. contains only pure + projections, but note that in this general situation the renormalized F{T) is no 
longer quadratic in T, as original F{T) was. 



2.4.1 Concrete Example 

To the clarify last statement, consider a simple example. Let V- pick up only terms of exactly the order e, 
^-^=(l7)e=o'andlet 

F{T) =T + beT^ + ce^T^ 



For such F(T), we can easily check that there is only one nonzero Q, namely Q{* * ) 
f = T- beT"^. Thus renormahzed F{T) is no longer quadratic in T: 

F{T) = F{f{T)) = T - 6eT2 + be{T - beT^f + ce^{T - beX'^f = 
= T + (c - 262)e2r3 + (63 _ ^cb)e^T'^ + Scb^e^'T^ - cb^e'T^ 



-be. It follows that 



(18) 



Naively this contradicts all the previous examples: with our previous experience we could rather expect F{T) — 
T + ce^T^. However, the point is that Q which would renormalize F in such a way, would contain 'P-positive 
terms and thus can not arise from the forest formula. The forest formula provides QiT) which satisfy the 
constraint (2). 

This is our first example, demonstrating that renormalized F can have very different shape - more sophis- 
ticated or simplified - from the original one. This implies that the second- level forest formula (5), adequately 
describing this shape modification, is not quite trivial. Let us show how it reproduces eq.(18). 

For ■F( « « ) we will have the same formula and result as in (15). For F( » » » ) we have (compare with (16)): 



Fr{- 



+2 I F(~)_l-(-P_)F(— ) + F{-^)—^{~V-) 



F{-^) 



F(^) 



F{~) 



1 



^) = F{~ 

i-V-)F{~) 



■) + 



= ce'^ + 2 {be{-V_)be + (-P_) [be{-V_)be] } 



{—V-) [be{—V-)be\ = and we obtain ■Fr( « « » ) = ce^ — 2b'^e^. In a similar way we can reproduce coefficients 

in front of T^, T^, . The interesting point is to understand, why FniT) vanishes for chain graphs with more 
than 6 vertices, as predicted by (18). In fact this is also almost obvious: it is easy to observe that trees with 
embedded box-graphs are not allowed, because V- in (5) acts upwards along the branches of the tree, and 
P_e" = for n > 1. All our boxes can contain only 2 vertices and the only nonzero F{T) arc those for chain 
graphs with 1, 2 or 3 vertices. Multiplication of 2 by 3 gives 6 - and for all the chain graphs with more than 
six vertices Fr^T) vanishes, in nice agreement with (18). 
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2.4.2 Another example 

To further deepen our understanding of problems which are solved with the help of forest formula, consider one 
more example with the same peculiar projector as in previous section 2.4.1 but a different 

F{T) = -T + beT'^ 

e 

In this case there is a single non- vanishing counter-term Q{T) : 

= -V. = -V-{be} = -be 

g(— ) = -be^ 

We see that V+Q{ * * ) = —be^ 7^ 0, what contradicts the condition (2). This is an example, where appropriate 
modification of (2), mentioned in footnote 2 is unavoidable: it is F{—— )Q(——-) rather than Q{——-) itself that 
should be "P-negative. This is indeed the case: 'P+F{-»—)Q{ « « ) = —V+be = 0. 

2.5 Loop diagrams, a need for vacuum energy: F{T) F{v, T) —v-\-T-\- eT^ 

Let us now represent T and by two loop diagrams and However, for both these diagrams F/F = • 

is a vertex of valence zero (no external legs), which is not represented in our function F{T). In order to make 
formula (4) applicable in this situation, we need to introduce an additional term into F, associated with this 
diagram. If we call the corresponding coupling constant (vacuum energy) v, then our new 

F(t;,T) = ?; + T + eT2 = . + O +e-0 (19) 

Now (4) can be applied, and the only non- vanishing Q(F) is 

q{0)=-V-{f[0)} = -^ (20) 

Since corrections to the T and v correspond to the graphs with two and zero external legs respectively, eq.(20) 
provides only a shift of v: 



v = v + Q[^^j =v-eT'^, f = T (21) 
so that renormalized F{v, T) =v-\-T + eT^ = t; -|- T is indeed independent of e. 

2.6 Mixed chain-loop diagrams, restoration of tensorial structures 

Let us now represent T by the loop O > but leave represented by the chain diagram « « . This time we 
have two bare couplings and 

F{v, T) = av + {a + b)T + ecT'^ = a» +a-<- +b Q + ec -»-<- (22) 

Now we have from (4): 

aQ (~) = -V- {f (— ) } = -ec (23) 

.q(o)^-p-{.(o)q(-)}-K-3-? (-) 

aQ {-^~) = -V- {2Q (— ) F (— ) } = 2^! (25) 
aQ{0)=-V-{F{0) Q(~)} = - V_ {Q(~)} = -2M!A (26) 
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To understand why the naively expected coefficients 2, 3 the numbers of ways we can cut the corresponding 
chain graph from the loop one, do not appear in (24) and (26), the tensor structure of v and T vertices should 
be restored. After this is done, our problem acquires the following form: 

av + af'i + 6f " + ecf'^f^i = av + aT'^ + feT" 

Due to the presence of two different tensorial structures this is actually a system of two equations: 

av + bf" = av 

Thus av = av — bT'^'-, and the loop graphs , , . . . are just the corresponding chain graphs « « , 



with ends glued together. 



For renormalization of v and T we obtain: 

aa aa 



+ 



V / V / aa aa 

q(„)„ + q(^„)^„+... = T-^T2 + 2(|)'t3-... (27) 



Since the values of Q[n(loop)j 

for the loop diagrams are determined by the values of Q[n\ for the corresponding 
chain diagrams, we can perform summation with the help of a generating function, in the same way as it was 
done in (11): 

F{v, f)=av + {a + b)f + ecf^ = a (^v + + . . }j + {a + b) (t - + . . . j + ec (t - ^T^ + ...)' = 

= av + {a + b)T 

Note, that even if a was originally put to zero, it acquires a-independent correction in renormalization process. 
The zero-valence vertex is automatically generated by the forest formula. However this happens in a singular 
way: the a-dependent correction is first divided and then multiplied by a — > 0. 

3 Generic polynomial from the valence-two vertex 

Now we are ready to attack a slightly more general problem than in the previous sections: take for a function 
F{T) a sum of two general polynomials with the same projector as in (7): 

F{T) = R{T) + eG{T) 
Renormalization problem reduces to the following equation: 

R{f)+eG{f)=R{T), 
and solution is represented as power series in e: 

^ = ^ ~ 'W) ^ - R"{T)G{T)) + ... (28) 

Immediate question is what is the origin of R' in denominator from the point of view of the forest formula. 
3.1 The simplest nontrivial F{T) ^T + aT'^ + eT^ 

We begin from example with the simplest nontrivial R: R{T) = T + aT"^ and G{T) = cT^. Substituting this 
into (28), we get: 

'l + 2aT+' (l + 2aT)3 ^^^> 
and expanding denominators in aT, obtain: 

f = T- e{T^ - 2aT^ + Aa^T^ - 8a^T*^ + ...)+ e^{^T^ - liaT^ + 42a^T^ + ...) + ... (30) 
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Now we attack this problem with the forest formula: 

F{T) = T + aT^ + eT^ = — + a — - + e ■ 

Q[3] = 0(— ) = -e, 



Q[4] = -7'_(2F[2]Q[3]) = 2ae, 
Q[5] = -V-{2F[2]Q[4] + 3F[3]0[3]) = 3e^ - 4a^e, 
Q[6] = -P-{F[2]Q^[3] + 2i^[2]g[5] + 3i^[3]g[4]) = 8ah - 13ae^ 

Thus for T we have: 

f = T-€T^ + 2aer^ + (Be^ - Aa'^e)T^ + (Sa^e - 13ae^)T^ + ...= 

= T- e{T^ - 2aT* + Aa^T^ - 8a^T^ + ...) + e^i^T^ - ISaT^ + ...) + ••• 

This is exactly the same as (30), but expanded first in aT and then in e. 

It is instructive to sum explicitly all the terms, contributing to the first and second order in e. As usual, 
this can be done with the help of generating functions Q{i){t) — J2n=3 ^^^Q{i)[M Q[2){t) = X]n=3 ^"<3(2) ['T']) 
where Q(i)[n] and Q(2)H denote the contributions to in Q[n] of orders e and respectively. 

Since Q[n\ does not contain terms of the order e°, only the terms with a single box-subgraph contribute to 
Q(i)[n + 1]: 

g(i)[n + l] = -P_(2i^[2]Q(i)[n]) 

For Q(i){t) this yields 

+3/-) [-51 f3 
0(i)(t) = Q[3]t' - 2atQ(,){t), Q(,){t) = 3^ = -e^^ 

For (5(2)[»^+l] we should take three contributions into account: F[2](5(2) Wi -P'[3]<3(i)M and i^[2](3(i)[n](5(i)[fc]: 
Q(2)[n + 3] = -V- ( 2i^[2]Q(2)[n + 2] + 3F[3]Q(i)[n + 1] + F[2] Q(i)[fc]Q(i)[n + 3 - fc] J 

V /c=3 / 

Therefore Q{2){t) satisfies 

Q(2) = - (2<F[2]g(2) + 3t2F[3]Q(i) + F[2]g^i)) 

and 

3et^g(^) + Qg(^)^ _ eH^ (3 + bat) 

^('^ " l + 2at ~ (1 + 2atf 

Thus, after re-summation, we reproduce (29). 

3.2 General F{T) 

To understand the way, in which the forest formula can differentiate polynomials and provide R'{T) in denom- 
inator in (28)), it is sufficient to sum up the terms of the order e^. Exactly as in previous subsection, this is by 
solving an equation for the generating function 



n=2 \n=l / 



Here R and G denote degrees of the polynomials R{T) and G{T) respectively, and for the sake of simplicity n 
was put to one. If 

Q(i)W = E*"Q(i)W' 

n=l 

then for n> R: 

Q(i)N = -'P-^ {2F[2]g(i)[n - 1] + 3i^[3]g(i)[n - 2] + 4^^[4]g(i)[n - 3] + . . . + RF[R]Q^i)[n -R+1]} 
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where is a projector which pics up only the e^-terms. The coefficients 2, 3, . . .R count the numbers of ways 
that we have to cut a box subgraph of corresponding the length from the chain graph [n]. When expressed in 
terms of Q(i){t), this recursion relation becomes: 



and finally 



g(i)(t) = - {2r2t + Brst^ + ... + RrRt"^-^) Q^^^it) - {g^t - g2t^ - ... - gct^) , 



R'{t) 



In a similar way we can perform term-by-term re-summations for all orders in e. 



4 Cubic equation. Valence-3 vertex and loops 

Let now T represent a valence-3 vertex. Then wc can use the forest formula to solve cubic equation T+eT^ = T, 
assuming that the projector is again given by (7). Solution is 

f = T-eT^ + 3e^T^ - 12e^T^ + 55e'^T^ + ...=T + ^(-)"c„e"T2"+i 

n=l 

where c„ are the numbers of rooted 3-valent trees with exactly n vertices, e.g. 



ci = 1 



C2 = 3 



+ 



+ 



(31) 



Wc now reproduce this solution from the forest formula. Like in the case of valence 2, there are different 
possibilities to represent our F{T) and thus different realizations of renormalization procedure. 

4.1 First realization 

The first way is to represent T + eT^ as 



F{T) =T + eT^= — «^ + ( 




(32) 



When we start calculating Q{T), the forest formula itself selects the relevant graphs from the huge set of all 
possible graphs with 4-valence vertices. 



Q 



-V- { F 








It is easy to notice, that since embedded graphs are not allowed, all non- vanishing Q(r) will have the form 

Q(r) = (-)^e^, 
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where k is the number of vertex in F. Thus the problem reduces to counting number of graphs in order k, 
which give nonzero contribution. However this becomes an easy task if we observe a one-to-one correspondence 
between 'triangle' graplis and rooted trees from (31): 



corresponds to 




When we grow three links on the some vertex of original tree, this corresponds to insertion of 
corresponding vertex of our graph, for example: 




into the 




Therefore we obtain from the forest formula the following result for T: 

( y\ { ( A ( A ( 



T = T + Q 





Q 





n=l 



4.2 Second realization 

The same function F{T) can be represented in a different way: 



F{T) = T + eT^ = 



+ 3e 



(33) 



This time we put additional coefficient 3, reflecting the symmetry of the problem. It becomes especially natural 
if we restore indices of a triple totally-symmetric vertex: T — > Tjjfe, so that (sums implied over repeated indices) 



J- ipq-^ mpq 



(34) 



and actually there are three different one-loop diagrams. 
The first counter-terms are: 



Q 



-V- 



— e 



(35) 



Note that there is exactly the same diagram at the r.h.s. and the l.h.s., no tripling occurs in this formula. 
At the next level there are diagrams of three different topologies: 



Q 




-W- { F 




Q 



3e^ 



(36) 



Q 



-2V- 



Q 



<<A<)-A-<A<)--A<MH-<A)]---'^ 



(37) 



(38) 



where coefficients 3 in (36) and 2 in (37) and (38) correspond to the number of possible choices of box-subgraphs. 
Altogether, taking into account three different orientations of the diagrams, we get for the counter-term 

Q{T): 

f = T- 3eT^ -|-3-(3 + 2-|-2x 2)e^T^ + ... 
and ^ 

F{f) = f + 3eT3 ={t- 3eT3 -|- 21e'^T^ + ...)+ 3e(^T - 2,eT^ + 27e^T^ -h • • • ) =T 
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5 Cubic equation. Two vertices: of valence 2 and 3 



We already encountered a problem with more than one coupling constant, but this was a rather trivial example. 
Now we reveal some more peculiarities of the forest formula's work with several couplings. Consider 



F{T) =T2 + e2Ti + T3 + esTi + 3eiTiT2 



+ €2 



+ 



+ £3 




(39) 



and V- pics up all the terms £3 £3 with all \a\ + > 1. 

Again, like section 2.6, due to the presence of different tensorial structures, eq. F{T) = T2 + T3 implies 
renornialization of two distinct functions: 

T2 + esTf = T2 
fs + eafi + 3e^Tf T2 = T3 

However this time these two equations are interrelated, and iteration procedure should determine T2 and T3 
simultaneously. Let us try to solve this system with the help of the forest formula. The first graph is simple: 



K>) = -''-{K^ 



-£2 



(40) 




> = ? 



However, when we look at the next graph with two external legs, we encounter a difficulty: 

/ A) 

«(-<I^)=-^- 



It turns out that we should first calculate 

/ A 

Q 

and only after that return to 




-V- i 




> = -£3 



(41) 



-V- I 




2e2e3 



(42) 
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The next graphs are: 




where (x3), (x4), . . . denote the numbers of possible positions of the dot or, in the case of (44), the number of 
orientations of the corresponding graph. Again note, that for (45) we need (40), for (46) and (47) we need (43), 
although these arc graphs with different number of external legs. 

The forest formula selects by itself the graphs which contribute to renormalization. It is because of connection 
between renormalization processes for T2 and T3 that we cannot determine first Q for graphs with two external 
legs, and then for those with three external legs - this should be done simultaneously, order by order: 

(40) (42) (46) (47) 

f2 = T2- £2 r| +2e2e3 Tl +^^2^ T2TI +2e2el T2TI +... (48) 

(41) (43) (44) (45) 

T3 = T3 - e3 T| -3e2 T2TI +'6el T| +3e2ei T| + . . . (49) 

To make clear the origin of every term, we keep reference to the graph, which it originates from. Now it is a 
simple exercise, to check that (48) & (49) solve original problem. 

Part II 

QFT examples 

QFT is a rather natural place to apply forest formulas - and it is what they were originally designed for. 
Perturbation theories in QFT often suffer from divergences, originating from momentum integrations in loop 
diagrams, i.e. from sums over tensorial indices in our diagram expansions. To handle divergent integral one 
applies regularizations, which introduce dependence on additional parameters, like an ultraviolet cut-off A in 
naive schemes, or e = 4 — d in dimensional regularization, or masses Mi in the Pauli-Villars approach. At the 
end of the calculation these parameters need be eliminated. This can not be achieved by simply throwing them 
away from the answer, because this would change mutual relations between different correlators, i.e. the shape 
of partition function. The risk is that after an arbitrary change we can loose possibility to obtain this partition 
function and these correlators from a single functional integral: they can fail to represent any kind of quantum 
mechanical system. Such risks are eliminated if one does not change the shape of partition function, but modifies 
its arguments, i.e. the bare couplings, which are exactly our parameters T from the previous sections. Forest 
formula gives us an expression for these bare couplings T in terms of the physical constants T: 

f = T + Q{T), 
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where Q{T) are singular counter-terms, which make partition function finite. In this sense the forest formula 
provides the resolution of rcnormalization pro})l{un in QFT. 

Forest formula operation does not depend on the nature of the function F(T) which we want to rcnormalize. 
The only thing one should care about to perform a self-consistent rcnormalization is the validity of the vertex 
criterium:. for every graph F with non-vanishing Q{T) there should be a vertex F/F in F{T), i.e. -F(F/F) 7^ if 
(5(F) 7^ 0. If this is true, rcnormalization in QFT is performed at the level of partition function, i.e. for all the 
correlators at once. However it is often possible to restrict consideration to correlators with particular external 
structure. 



6 0^ in one and two loops 

The simplest non-trivial example is provided by rcnormalization of the 4-point function in the cf>'^ theory at two 
loops (i.e. in the order T^). This example is still rather simple, nevertheless, it dmonstrates the main features 
and peculiarities of the forest formula. 
Our theory is: 

m = ^id,^r-^m^^^-^cl^' (50) 

We assume, that rcnormalization of field and of the mass term has already been performed, and thus we do 
not consider rcnormalization of kinetic and mass terms: valence-2 vertices and divergent graphs with 2 external 
lines are not considered (sec, however, sec. 9 in the Appendix below). We begin directly with the 4 point 
function, which is written as a sum of Feynman diagrams: 

Triple dots stand for omitted diagrams in the u and s channels. For Z{T\T) we take the following functions: 

Z (y<[ |t) = -iT, Z (1-loop diagramm|T) = {-iTf, Z (2-loop diagramm|T) = {-iTf, . . . 
Then the lowest coefficients ^'(F) are: 

^(x)=i. 

1 [ d^k i i 



(27r)4 F -m'^ {k+ pY - m 



Pauli-Villars regularization substitutes F^^'^^p^) by 

' 27 (27r)4 I^fc2-m2 fc2-M2y V(fc+p)2-m2 {k+p)^-M^J 



reg M\ 



327r2 



j-log-^+y^ dxlo, y^ + Oi-)\ = ^[-lo,^+I^,,i\p^)j (52) 



Similar formulas can be written for -F^^eo^M ^'^^ ^rea^\i ' them denote the squared momentum in the 



loop. Evaluation of F^'^'*^\p'^) = F ^^X^^^^^^^ straightforward, it is just the square of F^i'**): 
Another two-loop diagram is more complicated, moreover, it contains overlapping divergencies: 
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Here q is the momentum, running through the bigger loop, p is the sum of momenta in the bigger loop??? and 
Pi is the momentum of the right bottom external leg of the diagram. Substituting from (52) we obtain: 



p(2,t-(s+«)). ^ 



M2 



327r2 



+ 



[ ^ 

J (27r 



327r2 



Y -w? q'^ - J \ {q - pY - {q - p)^ - 

2 ( 



327r2 



(54) 



Since expression for every diagram naturally decomposes into two parts, one finite in the limit when M — > 
oo, and another, divergent, behaves like ^log^^ with some positive integer n, we ask V- to pick up the 



contributions like ( log 



with n > 0. 



Now we can apply the forest formula (4): 



Q 





= -V- 



F 



327r- 



■ log ■ 



M2 



(55) 



Q 



327r2 



= -V- { F 



2F 



Q 



(56) 



71^) - 7I^4:n^(^^) + 2 X ^ |- log ^ + /£)(p2)| log ^1 



log^ 



327r2 



(57) 




= -V- 




,3277^ 



l"g-^J -21og^-f/i„nP ) + ^/in )-2|-log-^+/)i;'(p )|log-^ 



2 / , ,r2 \ 2 



Collecting all these Q's into a single formula for T, we obtain: 

Such T indeed makes F{T) P-positive, and again, the resulting F{T) can be most straightforwardly obtained 
from the second-level forest formula (5): 

Fr{T) = F{f) 

(XX) - - (XX) - - (^) 4^ ^--^ - (XX) - (XX 
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327r2 







, , — x{l — x)p'^ 
dx log 5 



At the level of double-loop graphs we have to sum over the following forests for ^X0C3^ 



ra< xra :xH: 



and for 






According to the rules described in section 1.3, the forest /"(^^^CZX corresponds to 

{-V-)F 



(from now on we omit F =1)) while )(^3C3^ corresponds to 



{-v.){f 



{-V-)F 



If we represent the action of V- on the graph by surrounding it with an oval and denote the vertex which 
substitutes the contracted graph through ® , wc can symbolically write expression for the -R(r) = Ffi{r) as: 



R 



(59) 



«(XXX)=XXX 




(60) 



Now it is easy to recognize in i?(r) = -FR(r) the Bogolubov's R-operation [1, 2] and its graphical representation 
in eqs.(59) and (60). 

7 Renormalizable and non-renormalizable QFT example: in (i = 4 

and d > b 

In this section wc briefly illustrate the difference between renormalizable and non-rcnormalizablc QFT models. 
Note, that the forest formula works order-by-order in coupling constant, thus all graphs of a given order should 
be included. Coming back to eq.(51) for the cf)^ theory, we observe that in the third order there are more 
diagrams, for example. 



(61) 



It was omitted in eq.(51), because for d = 4 its Feynman amplitude does not contain divergent terms: 



< F 



= 



(62) 
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However, this statement depends on space-time dimension d, and this diagram diverges as Af* ^ for d> 6. For 
d = 5 this graph is finite, however, the diagram 

(63) 

docs diverge. This means, that in order to satisfy the vertex criterium and perform a self-consistent renormal- 
ization we have to include an elementary six- valence vertex Tg into F(T), in addition to the elementary 4- vertex 

After this the contributions from other graphs should be added, until all diagrams up to the third order will be 
taken into account. The diagrams to be added are, for example: 



Both of them diverge, since they contain a loop with two propagators, and they give contribution to Q{T). 
The first one contributes additionally to renormalization of Tg, however the second diagram generates Tg - an 
elementary eight-valence vertex: 

F(^)g(XX)=-^-{^(XX)}^° ^ 

In turn Tg generates Tio, and so on and so force. We see, that for d > 5, despite renormalization of F(T) is 
possible, an infinite number of vertices should be introduced to perform it self-consistently even in the third 
order in T. Such theories are called non-renormalizable. 

Situation in = 4 advantageously differs from higher dimensions. Due to (62), the diagram (61) does not 
contribute to renormalization of F(T) in (51) and generates no new vertices. And this property is not bounded 
to the order T^. It can be shown [l]-[3], that all divergencies in graphs with more than 4 external legs come 
from the sub-divergent 4-valence graphs. For example, although for d = 4 



this divergence does not force one to introduce the bare Tg vertex, because it is exactly canceled in the forest 
formula: 



F(*)«(^).P.{^(X 




0, 



and we can consistently keep F (^/"jC) =0. In the language of renormalized perturbation theory this means 

that all divergencies can be absorbed into a single coimter-tcrm for the 4- valence vertex: Q(r) which makes 
F[T) or any of the correlators in the theory finite, takes non-zero value only on the graphs with 4 external legs. 
Theories, possessing such property, are called renormalizable. 



Part III 

Appendices 

8 Proof of the forest formula 

After a number of concrete examples that were examined in the main part of this article, it should be a simple 
task to follow the proof of the forest formula. We remind that for a given function F{T) and linear projector 
V- it solves the equation (1), 

V-{f(t + Q{T)^} =0 (64) 

for the shift T ^ T = T + Q{T). In other words, it allows to modify the value of the function F in desired 
way, by shifting its argument while preserving its shape. This is important in physical applications, because one 
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can cure problems of partition function (e.g. divergences), without breaking its distinguished properties (say, 
keeping it inside the narrow class of integrable r- functions) . 

Representing 

F(T + Q(T)) = J2 FiT)Z{T\T + Q{T)) = ^ Fn{T)Z{T\T) 

r r 

we see, that since the functions Z{T\T) form a basis in the space of all functions of T, eq.(64) can be satisfied 
if and only if V-Fr{T) = for all T. (We remind that projector V- does not acts on the Z{T\T)). To find the 
coefficients Fji{r), we substitute Q{T) = Q{'y)Z{'y\T) and compare the two sides of equality: 

^F(r)z(r| T + ^g(7)z(7|T)) =^Ffi(r)z(r|r) 



Since 



Z(r| T + ^Q(7)Z(7|T)) =Z(7|T)+ ^ Q{j)Z{T\T)+ ^ Q(7i)Q(72)2(T|T) + . . . 

7 T,7:T/7=r T,7i,72:T/(7i-72)=r 



we get for Fr: 



Fr{T) = F{T)+ Y1 ^(r/7i-- -7^)^(71) •••0(7fc), 

{7iU...U7fc} 



(65) 



V-Fr{T) = 



V- < F(r) + F(r/r)Q(r) + ^' F(r/7i . . .7^)^(71) • • • Qin) \ = 0, 

{7iU...U7fe} J 

where prime in signals that we extracted the term with 7 = F from the sum. Now, supposing that 
V+F{T/T)Q{T) = 0, we arrive to the forest formula (4): 



F(r/r)g(r) = < 



F{T)+ ^^(r/71-- -7^)^(71) •••Q(7fc) 

{7lU...U7fc} 



Q.E.D. 

The second- level forest formula (5) can be obtained by solving recursion for Q{T) and substituting the result 
into (65). 



9 Kinetic terms, IPI diagrams and other peculiarities 

Here we briefly comment on some peculiarities, that were not considered in the somewhat oversimplified treat- 
ment of the section 6. 

The first point to explain is the situation with the IPI (one-particle-irreducible) graphs. They were omitted 
in the treatment of sec. 6 because they do not contribute to the counter-terms Q. This can be easily proved by 
induction. Indeed, consider a one-particle-reducible graph 




(external legs are not shown). Its Feynman amplitude is given by the product of propagator and two amplitudes 
for its IPI parts 71,72- Assuming by induction hypothesis that Q for all the one-particle-reducible graphs with 
smaller number of vertices is zero, we conclude that the set of box-subgraphs of F, which can potentially 
contribute to Q(F) is just the unification of the corresponding sets for its IPI parts, plus boxes which surround 
the whole 71 or 72. Thus wc get: 

Q(F) = V-{ (F(7i) - Q(7i)) (i^(72) - Q(72)) } X {— } = 
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because 

P_{f(7i) - Q(7i)} = 'P-{f{72) - 0(72)} = 

This derivation is clearly based on the fact, that all divergencies in QFT come from loops, while the one-particle- 
reducible graph is just the product of its components. 

The second point to clariiy concerns renormalization of mass and kinetic terms. To study this subject we 
need to introduce two more coupling constants in our treatment in 6 and, correspondingly, two more elementary 
vertices: 

the first one corresponds to the mass term and the second one to kinetic term. Lagrangian and elementary 
Feynman vertices look as follows: 

F{^) = h F{^)=p^, f(x)=1 

Z{^\T)=m^ Z{-^\T)=K, Z(X\T)=-^T 

Usually one assumes that the physical value k = 1, but its bare value k can be different. 

Now we can write the sum of all correlators which potentially can contribute to Q, i.e. are potentially 
divergent: 

=^ +^ +x +^ +^ +-e^ + (XX + -) + (>(^ + -) + 

Kxxx--)-(x5--)-(5x -■■)-■■ 

Since we have two different elementary vertices, corresponding to two different couplings, a new external struc- 
ture should be introduced to separate contributions to Q into two parts: Q for graphs with 2 external legs 
acquires additional label, equal to or 1. We have two projectors Pi"-* and T^i^-* for the self-energy graphs: 

V-^ picks up the singular terms, which are proportional to while picks up terms, proportional to p^. 
Accordingly, 

QW(r) = -pW{F(r) + ...} 

Now renormalization process may be done "as usual", 

self-energy 
graphs 

K=K+ 

self-energy 

graphs 

-if=-iT+ J2 

graphs with 4 
external legs 

For more elaborated treatment of renormalization with mass terms, see refs.[6, 8] for <p^ theory, ref.[5] for (p^ 
theory and refs.[8, 9] for QED. 
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